Orbital physics of polar Fermi molecules 
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We study a system of polar dipolar fermions in a two-dimensional optical lattice and show that 
multi-band Fermi-Hubbard model is necessary to discuss such system. By taking into account both 
on-site, and long-range interactions between different bands, as well as occupation-dependent inter- 
and intra-band tunneling, we predict appearance of novel phases in the strongly-interacting limit. 
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Creation of ultracold hetero-nuclear molecules opens 
the path towards experimental realization of strongly- 
interacting dipolar many-body systems. Depending on 
the constituent atoms, in moderate electric field these 
molecules can have large dipole moment of 1 Debye 
in their vibrational ground states [J-Q. In particular, 
fermionic molecules in presence of an optical lattice can 
be used to simulate various quantum phases, such as 
quantum magnetism and phases oft— J like models 
various charge density wave orders 0, bond-order 
solids ^ etc. One should also stress that in the strongly 
correlated regime, both in bosonic and fermionic systems 
the standard descriptions of single-band Hubbard model 
ceases to be valid. The effect of higher bands become 
important leading to novel phases like pair-superfluidity 

IMi. 



etc 



While most of the works have dealing with higher 
bands concentrated on bosonic systems, in this paper, 
we study dipolar fermions confined in 2D optical lattice 
Matt = Vo [siii^iTTx/a) +sm^{TTy/a)] + '-^z^, where Vq 
is the lattice depth, a is the lattice constant, m is the 
mass of the molecule, and is the frequency of har- 
monic potential in z direction. The dipoles are polar- 
ized along the direction of harmonic trapping. Usually, 
at low temperature and for low tunneling, the phase di- 
agram consists of different crystal states whose structure 
depends on the filling n [T]. In this paper, we derive a 
Fermi-Hubbard model for dipolar fermions including the 
effects of higher bands. We show that, even for moderate 
dipolar strength, it is necessary to take into account the 
excitations along the z direction. Simultaneously, in this 
regime, the interaction induced hoppings along the lat- 
tice give also important contributions. This changes the 
phases expected for a spinless Hubbard model including 
only a single band. Near n > 1/4, we find a crossover 
to one dimensional (ID) chains of decoupled Luttinger 
liquids along with density wave order. Near n > 1/2, 
we find that the system can be mapped to an extended 
pseudo-spin 1/2 Hubbard model with different emergent 
lattice configuration. We find a regime where chiral p- 
wave superconductivity emerges through Kohn-Luttinger 
mechanism with transition temperature Tc of the order 



of tunneling. This gives rise to an exotic supersolid, with 
the diagonal long-range order provided by the checker- 
board pattern of the lower orbital fermions, while the 
superfluidity originating from the fermions in the higher 
band. 

The Hamiltonian for the dipolar fermions in the sec- 
ond quantized form reads H — J d'^r5'^(r)iJo^(r) + 
i/d3rd3r'^'t(r)^'t(r')V(r- r')*(r')*(r), where ^'(r) 
is a spinless fermion field operator. In the units of recoil 
energy Eji — 7r^/i^/(2ma^), the single particle Hamilto- 
nian becomes Hq = — V^-|- Viatt(r)/£^fl and the long-rage 
interaction potential V(r) = D (l/r^ — 3z^ /r^) , where 
D = 2Trmd?/{lT?a) is a dimensionless dipolar strength, 
related to the electric dipolar moment d. For KRb 
molecules with a dipole moment of 0.5 Debye confined in 
the optical lattice with a = 345nm [l6[ one gets D ~ 8.6. 
We decompose the field operator in the basis of Wannier 
functions in the x, y directions and of harmonic oscillator 
eigenstates in z direction. For convenience we introduce 
orbital index a — {pml} denoting p, m and I excitations 
in X, y, and z direction respectively. In this basis the field 
operator *(r) = Y.i.a^i'y^i-yi.'^)^ where Wia{r) is the 
single-particle wave-function in orbital a localized on site 
i — ixBx +'iySy {^x and By are unit vectors in the proper 
directions). Fermionic operator ai„ annihilates particle 
in this state. The Hamiltonian can be rewritten in the 
following Hubbard-like form H — "hS^' -|- 
where 
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Parameters E^j and comes from the single parti- 
cle Hamiltonian and denote single-particle energy and 
nearest-neighbour tunneling in orbital a respectively. 
The inter-particle interaction has three contributions to 
the Hamiltonian (|lb[) : (i) the on-site interaction energy 
of fermions occupying different orbitals a and a' of the 
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FIG. 1. Parameters of the Hamitlonian for Vo = 8Er as 
functions of the lattice confinement, (a) On-site interaction 
— -jj^ (blue solid line -1-), and nearest-neighbor long-range 
interactions ^"qJ"'^ (solid black line -2-), (dashed 
black line -3-), and ^"'j^jf""'' (dotted black line -4-). The red 
line -5- shows the ratio Vp^^p^{2ex)/T^f^. (b) Magnitudes of 
the induced tunneling terms '^J — — (black solid line 

'^"''dj'^^^ (^^'-^ solid line -2-), and '^'"'"j^}^'^^^ (blue solid line 
-3-). The dashed black line -4- denotes the ratio tJ^^/Js for 
D = 10. 



same site Ug-a-i, (ii) the long-range interaction energy of 
fermions occupying orbitals a and a' of different sites 
Kto-'(^ ~ 3)1 (i'^'i) s-nd the tunneling from orbital a' at 
site j to the orbital a at site i induced by presence of 
an additional fermion at site i in orbital a" denoted by 

The Hamiltonian ^ is very general. To get a phys- 
ical understanding of its properties, we start by ex- 
amining the properties of density-density interactions. 
We calculate the interactions between few lowest bands: 
s = {000}, = {100}, py = {010}, p, = {001}, 
Pxz = {101}, and Pyz = {Oil}. We find that the on- 
site interactions Us^p^ = Us^py is always repulsive. It 
means that putting two fermions in s and Px or py band 
simultaneously is energetically unfavorable. Remarkably 
we find that /7s, is always negative. This surprising 
attraction stems from the presence of the fermionic ex- 
change term and the shape of the dipolar interactions (see 
supplementary material). Moreover, as it is seen from 
Fig. [T^, this interaction can not be neglected even for 
stronger confinements in z directions. For higher orbitals 
we find that Us.p, <Si Us,p^^ — Us,py^ < 0. In addition for 
long-range interactions we find that Vs^s{i) > Vs,p^{i) > 
^2,Pz(*) > 0- From this analysis we conclude that for 
polar molecules there always exists some critical dipo- 
lar strength for which single-band approximation breaks 
down since two particles can occupy the same site. This 
critical behavior is controlled by the on-site energy cost 
A = Ep^ + Us,p^- In this paper we consider situations 
when A is positive. 

Next we discuss the role of the interaction induced 
tunnelings in Hamiltonian ([I}. Counter intuitively the 



most important contribution does not come from the 
induced tunneling in pz band {Tp^^p^), but from the 
inter-band tunneling which changes Pz orbital to the 
Pxz and Pyz ones (T^^ (e^,)). Note, that this inter- 
band tunneling is absent for usual single-particle tunnel- 
ing due to the properties of single particle Hamiltonian 
[l9| . From properties of p orbital states it follows that 
Tp^ p^^{-ex) = -Tp^^p^^{ex)- The relation of this term 
to other interaction-induced tunnelings is shown in Fig. 
[TJa. From the above analysis we introduce simplified, but 
realistic model of polar Fermi molecules confined in 2D 
optical lattice by taking into account effects of interac- 
tions between orbitals tr £ {s,Pz,Pxz,Pyz}- 

Since we consider the case when A > 0, therefore 
for low filling all fermions occupy only the s orbital 
states. For filling 71 = 1/4 and large enough D (> 3) 
there is non vanishing single-particle excitation gap and 
the system is in the s-band insulator state (denoted by 
blue spheres in Fig. [2^) Q- Situation change dramat- 
ically for higher fillings. It can be simply understood 
with energy arguments. The energy cost of putting ad- 
ditional particle in the vacant site is given by i?vac = 
Vs.s{ex) + 2Vs,s{ex + By) -I- . . .. In contrast the cost of 
putting additional particle to the Pz orbital of an occu- 
pied site Eocc = A -I- 2Vs.p^ {^Bx) + . . .. For D larger than 
some critical strength one finds that i?occ < ^'vac- As 
an example, such conditions are fulfilled for Vq = SEn, 
D — 10, and hH. < 14i?fl. Consequently, additional par- 
ticles start to fill Pz band of previously occupied sites. 
In this scenario energy conserving dynamics of the sys- 
tem comes from the second-order processes involving tun- 
neling to the next occupied site (along x direction in 
Fig. [2^). To the leading order, this effective tunneling 
is given by 

tL « KpAex f/m.^pA + Ep^^ - EpJ. (2) 

Thus, the pz fermions will only move along one direc- 
tion chosen by the insulator checkerboard geometry in 
s-band. The resulting system can be thought as stacks of 
independent, one-dimensional chains without inter-chain 
tunnelings. Each chain is characterized by the effec- 
tive tunneling tJ^ j and the nearest-neighbor interaction 
Vp^^p^{2ex). Relation between these parameters can be 
tuned experimentally by changing the trap confinement 
(red line -5- in Fig.[T^) as well as the lattice depth. More- 
over, due to the absence of the inter-chain tunnelings, we 
can apply Jordan- Wigner transformation to each chain 
to map it to chain of independent spin 1/2 XXZ model 
resulting in a crossover to one-dimcnsional system. For 
22^iff > ^p^.pA'^^x) such a model will show metallic be- 
haviour with the low energy sector is dominated by the 
bosonic excitations. From Fig. [T] we see that indeed 
V^_p^ (2ej;)/rjjf < 2 for any confinement for Vq — 8Eji 
and _D = 10. Thus we will see metallic behaviour of the 
system due pz fermions with the s orbital fermions form- 
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FIG. 2. Pictorial diagram of the difTerent checker board 
lattices. The blue and shaded-red sphere denotes s-orbital 
fermions and the smaller black sphere denotes pz orbital 
fermions. (a) Checkerboard lattice at n = 1/4 filling. The 
Pz fermions will move with effective tunneling tJ^j (arrow 1) 
only along the shaded region making a stack of ID chains. 
Interaction between neighboring p-band fermions is equal to 
T^^,p^(2ea;) (arrow 3). (b) Checkerboard lattice at n = 1/2. 
The blue and thick-red line constitutes two different sub- 
lattices which are not connected via tunneling processes since 
the tunneling T^^j (arrow 2) is much smaller than t\^. (c) 
Density- wave structure at filling &n = \ /4. The dark and light 
shadings denote higher and lower density of the pz fermions 
respectively. 



ing a checkerboard density pattern as shown in Fig. [5^. 

Next let us discuss the case of filling n = 1/2. Then 
for low dipolar strength D, due to the same reasons as 
before fermions will occupy only the s-band. To look into 
properties of the ground state, we express the average 
density (ni,s) = (1 + (— 1)*"=^*"(5)/2, where 8 is the order 
parameter. We also define the single-particle imaginary 
time Green functions Q{% — J,t) = (Ta^ s(T)at ^(0)). By 
following the procedure described in [13, [l3 we find the 
following equations for Q in the momentum space 



[cj + ^ - 2V{\ - 5')] 01(^,0;) - efcG'2(fc,w) = 1, (3a) 
[cj + ^-2y(l + (5')]G'2(fe,w)-efcG'i(fe,w) =0, (3b) 

where in the position space G\{i) {G2{i)) is equal 



to G{i) for 



even (odd). The kinetic en- 



ergy eic = 2 Js {cos kxtt + cos kytt), the effective poten- 



tial V 



and S' = 5(E, 



Sodd 



J2o=£ieeven^ssi'>')) /V . Thcsc mean field equations for Gi 



and G2 are similar to the ones found for extended Hub- 
bard model, with a renormalized nearest neighbour inter- 
action and density imbalance IT], ll8| . In our case we con- 
sider the neighbors within a radius of |i| < 5. Then, by 
solving equations (|3]), we find that in the strong coupling 



limit (5=1 — 



3J; 



{d ^ .98 for previous parameters). 



Thus, we find that in this limit the ground state of the 
system is the checkerboard insulator (blue and thick-red 
spheres in Fig. 

Next we look to properties of the system with ad- 
ditional particles. To do it we define deviation from 
half- filling Sn — n — 1/2 and we introduce correspond- 
ing chemical potential fj,{dn). From energy arguments 



we find that two scenario can happen. The additional 
fermion (i) occupies a vacant site with energy cost 



4K. 



8I4,s(2ea; -I- By 



or (ii) it goes 



to the pz orbital of an occupied site with energy cost 
Eocc = A -I- Ws^p^ie^ + By) + T4,p,(2ej;) + Conse- 
quently, in the second scenario (when Eocc < -E'vac), all 
extra fermions will occupy the Pz orbitals of the already 
occupied sites. As an example, such conditions are ful- 
filled for Vo = 8Er, D = 8, and hft < WEr. In such 
a case Sn corresponds to the filling of Pz band fermions. 
The parallel tunneling of the Pz fermions between the 
occupied sites will again arise from the second order pro- 
cesses Moreover, tunneling to the diagonally occu- 
pied site T^if --[Js- T;^,pMx)]^ /\Us,pJ for D ^ 10 
it is 400 smaller than T^ff. Consequently fermions in 
the Pz orbitals can move in independent square sub- 
lattices (either the thick-red or blue sub-lattice shown 
in the Fig. [21d). Note, that fermions can not tunnel be- 
tween different sub-lattices. Thus we can describe the 
system of the pz fermions in the blue (thick-red) lattice 
as pseudo-spin up (down). By introducing operators Cjs, 
where s E {t, 1} the effective Hamiltonian can be written 



as -ffe 



iff '^{ij} ^is'^js 



Hi„t. with 



« {ij} 



(4) 



where fiis — cl^Cis- For convenience we introduce Vf-f = 
Vp^^p^{2ex) and Vfi — Vp^^p^{ex + ey). Note, that now 
{.} is understood as a nearest- neighbor in a given sub- 
lattice. Nearest-neighbors between different sub-lattices 
is denoted by [.]. The modified lattice constant of the 
sub-lattices is a = 2a. In this way we are able to investi- 
gate the system properties with weak-coupling theory. To 
do this we transform to momentum space and introduce 
charge fiuctuations pq — ^ '^'k+q s'^k,s and spin fluc- 
tuations Sq = g sc^^g ^Ck.s operators. Subsequently 
we write 

^^int ^jYl C^^nVq + VuPq) PqP~q 



+ \T.(^^nVq~Vul3q)SqS^q, 



(5) 



where r]q = 2{cos{qxa) + cos{qya)) and /3q = 
4(cos((7a;a/2) cos((7j,a/2)). The system with interactions 
described by ([5]) can manifest three possible magnetic 
instabilities: charge-density wave (CDW), spin-density 
wave (SDW) , and ferromagnetic instability. At Sn = 1/ A 
(each sub- lattice is half- filled in p orbital), as /3q = 
for nesting vector qa = (±7r, ±7r), interaction in the spin 
channel becomes repulsive and therefore SDW order is 
absent. We also find that Stoner ferromagnetism occurs 
only for very low temperature on the order of 10~^ x rj^^ 
and therefore it can hardly be observed. 
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Finally, we examine the CDW instability at Sn = 1/4 
with modulation qa = (±7r,±7r). Within weak coupling 
theory, the transition temperature of the CDW state for 
the system described by the Hamiltonian Q is given 



(a) 



(b) 



by Tci 



ST'lfCxp 



T^tf/Vf^ )• It means that 



transition temperature is of the same order as tunnel- 
ing energy in the s-band. For example, when Vq = 8i?ij, 
hfl = IOEr, and D = 8 we find that Tcdw ~ 0.35J^ 
(~ 2 nK) . One should note that due to the relative shift 
of sub-lattices the resulting density modulation in this 
phase looks like stripes rather than standard checker- 
board structure as shown in Fig. [2j:. 

Lastly, we investigate the emergence of triplet su- 
perconductivity between the same pseudo-spin fermions, 
arising via Kohn-Luttinger mechanism (KL) [l^. We 
look for Cooper pairs with chiral p-wave symmetry. The 
effective interaction between fermions in KL mechanism 
in terms of the scattering momentum k — k' = q can be 
written as 



Kff sAq) = - E [(^tt 4 + f^l) Qi^p (6) 
p 



where Qq^p = ■• fi^) is the Fermi distri- 

bution function and ep — 2T^^{cos{qxa) + cos(g.ya)) 
is the dispersion. The summation in (IS3p comes 
from taking into account the second-order terms rep- 
resented by diagrams shown in Fig. [3^. By perform- 
ing the integration over the momentum in the limit 
of T — >■ 0, we finally get antisymmetric part of effec- 
tive couphng {Vett{q)}_ = —X{T, fi){sin{kxa) sin{k'Ji) + 

sin{kya) sm{k'ya)) where A(T,^) = 2V^^ + -^Fi{T, fi) - 
-F2{T,fi). Functions Fi and F2 originate in the 



second-order corrections and their detailed forms are 
given in the supplementary material. The point is that, 
due to the Van-Hove singularity in density of states, func- 
tion F2 contains a logarithmic divergence. At the same 
time function Fi is analytical due to the dressing of the 
density of states. This means that there always exists 
finite critical ^ above which the interaction is attractive 
and superfiuidity appears. From the BCS theory one can 
get an estimate of the transition temperature Tc (deriva- 
tion is shown in Supplementary material) . In Fig. ^jp we 
plot the transition temperature Tc as a function of de- 
viation Sn for example parameters discussed previously. 
For 5n - 0.22 we get - 0.2 1 nK). At half-filhng 
with 5n — 1/4 for the red and blue lattices, the transition 
temperatures for the CDW and superfluidity are similar. 
This will result in a competition between both instabil- 
ities. A detailed account of such scenario is beyond this 
scope of this paper. 

In conclusion, we have derived a generalized Hubbard 
model for dipolar fermions in an optical lattice by tak- 
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FIG. 3. (a) Diagrammatic representation of the second or- 
der contributions in (|S3[) . (b) p-wave superfiuid transition 
temperature Tc as a function of density n = 1/2 + Sn. 



ing into account higher orbitals. We have shown that 
the effect of these higher orbitals leads to new phenom- 
ena. Due to the strong interaction-dependent hopping 
terms in higher orbitals, these systems can be described 
by effective weakly-interacting theories. For particular 
parameters, near n > 1/4, we found a cross-over to the 
one-dimensional physics resulting in simultaneous metal- 
lic and density wave properties. For other set of param- 
eters, n > 1/2, the system can be described by a weakly 
interacting Hubbard model with pseudo-spin originating 
from the lattice geometry. Using the Kohn-Luttinger 
theory, we found a transition to the chiral p-wave su- 
perfluidity without destroying the checkerboard order. 
The parameters used here are currently experimentally 
achievable. 
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Supplementary Material for Orbital physics of po- 
lar Fermi molecules 

DERIVATION OF THE PARAMETERS U,,p^ , 

-^Pz.Pz -^Pz,Pxz 

In this section we represent the on-site interac- 
tion term Us^p^ and interaction-induced hopping terms 
'^Pz,Pz terms of the single-particle wave- 

function Wia{r) in orbital a localized on site i. Or- 
bital index a = {pml} denoting p, m and I levels 
in X, y, and z direction respectively. Then the s or- 
bitals can be written as Wis{r) — Wi^o{x)wiyo{y)4>o{z) 
where Wi^o{x),Wi^o[y) are the lowest band Wannier 
functions and (l>o{z) is the lowest orbital of the har- 
monic oscillator in the z direction. Similarly we can 
write Wip^{r) = Wi^o{x)wi^o{y)(j)i{z) and Wip^^{r) = 
■Wi^i{x)wi^o{y)(j)i{z). Here Wi^i{x) is the Wannier func- 
tions in the first band and (j)i{z) is the first excited state 
of the harmonic oscillator in the z direction. For simplic- 
ity we took j = i + Bx- From this we can write various 
parameters as, 



Us,p. = J {'WiAx)w^yo{y)V{'w^Ax>^yo{y')V^l,oiz. z')V{r - r')drdr' (Sla) 

TLpzi^-) = J w;,.o(a;)«;i.o(a;Xo(a;'){«'i„o(y)«^i„o(2/')}'*i,o(^,^')V(r - r')drdr' (Sib) 

TLp.M-) = /«a.i(-^)"'*.o(a:Xo(2^'){«^<.o(y)^^.„o(y')}'*i,o(^,^')V(r - r')drdr' (Sic) 

$i,o(z, z') = {(6i(z)}2{/(z')}' - <^i(z)/(z)0i(z')'^"(^')- (Sid) 

The integrations over z, z' can be done analytically using convolution theorem in the momentum space. Consequently 
in the momentum space we get 



V{k±) = '^hoiz, z')V{r - r')dzdz'^ 



[kdy - \l^kd{l + (fcrf)^)erfcx 



/ kd\ 



(S2) 



where denotes Fourier transform, k = \k±\ = 

^Jk-^c + ky, and Iz — {h/mQy^^ is a natural harmonic os- 
cillator length unit. It is important to note that V{k±) 
is always negative for any A; j_ . This explains the appear- 
ance of the attractive on-site interaction for any value 
of the confinement along the z direction. Also we define 
erfcx(a;) = exp(a;^)erfc(a;) where erfc(.) denotes comple- 
mentary error function. 



DERIVATION OF EFFECTIVE INTERACTION 
IN THE TRIPLET CHANNEL 



In this section we derive the p-wave interaction from 
the Kohn-Luttinger effective interaction in terms of the 
scattering momentum q = k — k'. For this purpose we 
rewrite Eq. (6) from the original paper 



(S3) 
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where Qp,q = e^p—q^l^p) with /(.) being the Fermi distribution function. First we put the expression of r]q = 
2{cos{qrcd)+cos{qya)) and f3g — 4(cos(ga;a/2) cos(gya/2)) back to (|S3p . As we are interested in p- wave interaction, after 
expanding (jS3p in terms of the momentum kx, k'^ ,ky, k'y, we keep terms proportional to sin fca;a sin fc^a+sin fcj,asin kyli. 
In this way we get 

{Kff (g)} = 1^2% - Yl 0(9.P)[4(%)' - 8(%)2{cos(fc,a - px~a) + cos{ky~a - Py~a)}Yj (sin fcajQ, sin fc^a + sin fcj,a sin k'yCi) 

+ 4(T^^)^ ^^(5(q,p)(sinfc2;asinA;^asin^Pa;a + sinA;j,asinfcyasin^Pya) (S4) 
p 



where q — k — k'. By converting sums to integrals in 
(jS4p we have to compute terms of the form J dpQ{q,p), 
J dpsin^ pxabQ{q,p), and J dp cos PxabQ{q,p). Func- 
tion Q{q,p) is a continues counterpart of Qq^p. In the 



limit of vanishing temperature T — > 
integrals / dpG(p)^^^^£^i^^^^^ta^ 



we approximate all 



e(p— q)~e(p) 

for arbitrary function G{p). The effective density of 
states reads Ncs{e,^p) — J dkG{k)S{e — €k). We see 
that the only first term inside the third bracket in (|S4p is 
not dressed by cos(pa:a) or sin(pxa). Hence the effective 
density of states for this term contains Van-Hove singu- 
larity. All other terms with in (IS4[) . due to the dressing 
by cos(p2;a) or SYii{pxd), are analytic. For convenience, 
we re-express {Kff (?)}_ = ^-^(r, ^)(sin(fc2,a) sin(fc^a) + 

sin(A:j,S) sin(fc;S)), where A(T, fi) = 2%-(--%Fi(T, fi)~ 
—^-F2{T,fj.). The functions Fi and F2 are given by 



where F{x) — E{x) — (1 — x)K{x) with K{.) being the 
Elliptic integral of the first kind and E{.) is the Elliptic 
integral of the second kind. Then we can write the BCS 
equation for the transition temperature Tc as 
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where the effective density of state is given 

^ - ek) SVC? [kxa-b) ~ 



by 7Ve//(e,M) = EkK^ 



eff- 



Fi = -F(l-(T+H)V(4Til,)2) (S5a) 

F2 = ^x([i-(r+H)V(4rii,f), (S5b) 
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